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Synchronization of spatiotemporally chaotic extended systems is considered in the context of coupled one- 
dimensional Complex Ginzburg-Landau equations (CGLE). A regime of coupled spatiotemporal intermittency 
(STI) is identified and described in terms of the space-time synchronized chaotic motion of localized structures. 
A quantitative measure of synchronization as a function of coupling parameter is given through distribution 
functions and information measures. The coupled STI regime is shown to dissapear into regular dynamics for 
situations of strong coupling, hence a description in terms of a single CGLE is not appropiate. 
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Two issues of high current interest in the general field of 
nonlinear dynamics are the quantitative characterization of 
different regimes of spatiotemporal complex dynamics in ex- 
tended systems [|l]] and the synchronization of chaotic oscilla- 
tors iQ]. The characterization of low dimensional chaos is now 
a mature subject with well established techniques, including 
techniques of chaos control. In this context, the demonstration 
that the familiar phenomenon of synchronization of two regu- 
lar oscillators |Q| by a weak coupling can also be displayed by 
chaotic oscillators is an important new idea. This conceptual 
development has opened a new avenue of research with inter- 
esting practical implications. Chaos in extended systems is a 
much less mature subject and many investigations are still at 
the level of classifying different types of behavior Concepts 
and methods of Statistical Mechanics are commonly invoked 
in terms of "phase diagrams" and transitions among different 
"phases" of behavior [Q-^. Still, the possibility of a synchro- 
nized behavior of spatially extended systems in a spatiotem- 
poral disordered phase is an appealing idea that we address 
in this Letter. More specifically we will consider an extended 
one-dimensional system in a chaotic regime known as Spa- 
tiotemporal Intermittency (STI) and we will characterize a 
coupled STI regime. 

By synchronization of two chaotic oscillators Oi and O2 it 
is meant in a strict sense that plotting the time series Oi(ti) 
vs 02{ti) one obtains a straight line of unit slope. For many 
practical applications, synchronization of chaotic oscillations 
calls for an expanded framework and the concept of "general- 
ized synchronization" has been introduced [||,^ as the appear- 
ance of a functional dependence between the two time series. 
In this context we understand here by synchronization the sit- 
uation in which Oi(ti) becomes a given known function of 
02{ti), while for independent chaotic oscillators Oi{ti) and 
02{ti) are independent variables. Transferring these concepts 
to spatially extended systems, we search for correlations be- 
tween the space(a;i)-time(fj) series of two variables Oi {xi , tj ) 
and 02{xi,tj). The synchronization of Oi and O2 occurs 
when these two space-time series become functionally depen- 



dent. This idea is different from the one much studied in the 
context of coupled map models in which the coupling and 
emerging correlations are among the local oscillators of which 
the spatially extended system is composed. Here we search 
for correlations of two variables at the same space-time point. 

Our study has been carried out in the context of Complex 
Ginzburg Landau Equations (CGLE) which give a prototype 
example of chaotic behavior in extended systems [|lO|jll]] . Our 
results show that the coupling between two complex ampli- 
tudes Ai and A2 (Oi = \Ai\ and O2 = |^2|), in a STI 
regime described below, establishes spatiotemporal coiTela- 
tions which preserve spatiotemporal chaos but lead to a syn- 
chronized behavior: Starting from the independent STI dy- 
namics of Ai and A2, coupling between them leads to a STI 
regime dominated by the synchronized chaotic motion of lo- 
calized structures in space and time for Ai and A2. An ad- 
ditional effect observed in our model is that the coupled STI 
regime is destroyed for coupling larger than a given threshold, 
so that the two variables remain strongly coiTelated but each 
of them shows regular dynamics. At this threshold maximal 
mutual information and anticoiTelation of \Ai\ and 1 are 
approached. 

The CGLE is the amplitude equation for a Hopf bifurcation 
for which the system starts to oscillate with frequency utc in a 
spatially homogenous mode. When, in addition, the Hopf bi- 
furcation breaks the spatial translation symmetry it identifies 
a preferred wavenumber Kc- In one-dimensional systems the 
amplitudes Ai and A2 of the two counterpropagating travel- 
ing waves with frequency ujc and wavenumbers ±A'c becom- 
ing unstable satisfy coupled CGLE of the form 



dtAi^2 = + (1 + ia)dlAi^2 

-(l + z/3)(|Ai,2p + 7|^2a|') Au 



(1) 



Eq. (|l]) is written here in the limit of negligible group velocity. 
In particular, this limit is of interest to describe the coupled 
motion of the two complex components of a Vector CGLE. In 
this context, (|l]) is used to describe vectorial transverse pat- 



1 



tern formation in nonlinear optical systems, and Ai 2 stand 
for the two independent circularly polarized components of 
a vectorial electric field amplitude [|l2|jlj]. The parameter /i 
measures the distance to threshold and 7 is the coupling pa- 
rameter, taken to be a real number. 

Homogeneous solutions of Eq. ([l|) are of the form 



(2) 



where: lji,2 = -P{Qi,2 + iQi.i)- For 7 = 0. Qi,2 
and the two amplitudes satisfy independent CGLE whose 
phase diagram has been studied in much detail in terms of 
the parameters a and f3 [p4|,[7||. For 7 = 0, solutions of the 
type (^, and other plane waves of different periodicities, are 
known to be linearly stable below the Benjamin-Feir (BF) line 
(1 + af] > 0). Above this line regimes of phase and defect 
chaos occur. However, for a range of parameters below the 
BF line there is an additional attractor, coexisting with the one 
of plane waves, in which the system displays a form of spa- 
tiotemporal chaos known as STI. In this attractor the solution 
is intermittent in space and time. Space-time plots of |yli| or 
1^2! in the STI regime for 7 = are qualitatively similar to 
the ones shown in Fig. |l](top). The question we address here 
is how the STI regimes of Ai and A2 change when the cou- 
pling 7 is introduced. We first recall that for a weak coupling 
situation (7 < 1) the solution (||) with Qi 2 = m/(1+7) is lin- 
early stable below the same BF-Hne 1 + a/3 > [|l|] whereas 
the solutions with Qi = or Q2 = are unstable. For large 
coupling 7 > 1 the competition between the two amplitudes 
is such that only one of them survives, so that linearly stable 
solutions are either Qi = ^/JI, Q2 = 0, or Q2 = ^/JI, Qi = 0. 
For the marginal coupling ^ = 1, Q\ + Q\ = n and the phase 
X = arctan(Qi/Q2) is arbitrary. In addition to these ordered 
states we also find a STI attractor for coupled CGLE and val- 
ues of a and /? which are in the STI region of a single CGLE. 
Changes of such STI behavior with varying 7 are shown in 
Fig. I?]©- 

For small coupling (7 << 1) we observe that \Ai \ and | A2I 
follow nearly independent dynamics with the flat grey regions 
in the space-time plot being laminar regions separated by lo- 
calized structures that appear, travel and annihilate. In the 
laminar regions configurations close to (H) with Qi ~ Q2 oc- 
cur. Disorder occurs via the contamination by localized struc- 
tures. These structures have a rather irregular behavior and, in 
a first approach they can be classified in hole-like or pulse- 
like JTl]]. In Fig. U these hole-Uke and pulse-like structures 
are associated to black and white localized structures respec- 
tively. It is argued that the domain of parameters in which STI 
exists in the limit 7 ^ is determined by the condition of 
stability for those localized structures [ 16 1. As 7 increases we 
observe two facts: First, both \ Ai \ and IA2I continue to dis- 
play STI dynamics although in larger and slower space-time 
scales. Second, and more interesting, is that the dynamics of 
1^1 1 and 1^2! become increasingly correlated. This is easily 
recognized by focusing in the locaUzed structures: A black 
traveling structure in the space-time plot of \Ai \ has its cor- 
responding white travehng structure in the space-time plot of 



1^2 1 and vice versa. This results in laminar states occurring in 
the same region of space-time for \Ai \ and |^2 1- The coupled 
STI dynamical regime is dominated by localized structures in 
which maxima of |^i | occur always together with minima of 
1^2 1 and vice versa. In the vicinity of the localized structures, 
and emerging from them, there appear travelling wave solu- 
tions of (|l]) but with different wavenumber for \Ai \ and 1^2 1 
so that I All 7^ 1^21- Eventually (going beyond the marginal 
coupling 7 = 1) the STI dynamics is destroyed and \ Ai \ and 
1^2 1 display only laminar regions, in which either \Ai \ or 1^2 1 
vanish, separated by domain walls. 




FIG. 1. Space-time plot ot the modulus \ A\ \ (left) and \A2 \ (right) 
for four values of the coupling parameter 7: From top to bottom, 
7 = 0.1,0.5,0.95 and 1.05. The horizontal axis represents space 
and the vertical axis time (2000 time units for 7 = 0.95, 100 for 
7 = 1.05 and 200 in the other two plots). The grey levels change 
linearly from the minimum (black) to the maximum (white) of the 
modulus. The parameters are /i = 1, q = 0.2 and 13 = —2.0. 
For 7 = 0.1, the structures appearing in \A\\ and \A2 \ are almost 
independent. However, for 7 = 0.5, the defects in (black and 
white structures) becomes slightly synchronized with the defects in 
1^2! and, for 7 = 0.9, the synchronization is almost complete. 

In the optical interpretation of ([l]) the laminar regions with 
1^1 1 = 1^2! corresponds to transverse domains of linearly 
polarized light, although with a random direction of linear po- 
larization. The localized structures are essentially circularly 
polarized objects since one of the two amplitudes dominates 
over the other. Around these structures the plane wave solu- 
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tions with \Ai\ ^ \A2 \ have different frequencies, so that they 
correspond to depolarized solutions of ([]) [p^. As 7 > 1, 
localized traveling structures dissapear and one is left with 
circularly polarized domains separated by polarization walls. 




FIG. 2. Comparison of the states reached at long times start- 
ing from noise for (from left to right, and from top to bottom) 
7 — 0.1,0.5,0.95 and 1.05. The joint probability distribution 
p(|^2|, 1^1 1) is shown as a 3D surface. The vertical scale is arbitrary, 
the same in the first three plots and three times larger for 7 = 1.05. 
On top of each surface, |Ai(a;,t)| vs |^2(a;,t)| are shown in the 
form of a dotted plot obtained from the values of |^i| and \A2\ at 
space-time points during a time interval of 50 units. 

It is usually argued that for 7 > 1 the dynamics of the 
coupled CGLE (|l]) is well represented by a single CGLE since 
only one of the two waves survives. This is certainly not true 
in the STI domain of parameters considered here since a single 
CGLE would give rise to STI dynamics whereas the coupled 
set (|l]) does not for 7 > 1. In general a description in terms 
of a single CGLE would not be reliable for parameter values 
at which the single amplitude dynamics produces amplitude 
values close to zero. 

We next show that the correlations observed for increasing 
7 in Fig. |l] are in fact a kind of spatiotemporal synchroniza- 
tion, in the generalized sense defined in ||8|^]. To this end a 
characterization of the synchronizing process can be given by 
analyzing the joint distribution of the two variables. This dis- 
tribution and values of \Ai\ versus |j42| are plotted in Fig.|[ 
The cloud of points correspond to the different space-time 
points of Fig. |l|. For 7 << 1 we obtain a diffuse cloud 
of points indicating essentially independent dynamics. The 
concentration of points around = 1^2 P = m/(1 + 7) 

corresponds to the laminar regions, but excursions away from 
that solution are independent. As the coupling is increased 
with 7 < 1 the cloud of points approaches the curve given by 



1^1 



1^2! = This indicates synchronization of the dy- 



namics of structures departing from the laminar regions. The 
points with larger values of \Ai\ and smaller values of \A2\ 
(and vice versa) correspond to the localized traveling struc- 
tures. Intermediate points among these ones and those around 



1^1 1 = 1^2! correspond to the regular solutions of non-zero 
wavenumber that surround the localized structures. The spe- 
cial case of marginal coupling is discussed below, but as we 
enter into the strong coupling situation (7 > 1) the cloud of 
points concentrates in the regions \Ai\'^ ~ fi, \A2\ = and 
1^2 P = 1^1 1 = 0, while intermediate points correspond to 
the domain walls separating these ordered regions. It should 
be pointed out that we are considering just the modulus of the 
complex fields ^i,2- The coupled phase dynamics does not 
show synchronization, at least not in an obvious manner, so 
that we are in a case of partial synchronization as considered 
in [|17]1. 

A quantitative measure of the synchronizing process can 
be given in terms of information measures [|l8|]. The en- 
tropy H{X) — — lnp(a;), where p{x) is the 
probability that X takes the value x, measures the ran- 
domness of a discrete random variable X. For two ran- 
dom discrete variables, X and Y , with a joint probabil- 
ity distribution p{x, y), the mutual information I{X, Y) = 

-J2x.yPi^^y)^^[pi^)piy)/pi^^y)] S.^^^^ ^ measure of the 
statistical dependence between both variables, the mutual in- 
formation being if and only if X and Y are independent. 
Considering the discretized values of \Ai \ and |yl2| at space- 
time points as random variables X = \Ai\, Y = \A2\, their 
mutual information is a measure of their synchronization. In 
Fig. |](left) we have plotted the mutual information and the 
entropy of |^i| and \ A2\ as a function of 7 Jl^]. This graph 
shows that the entropy of |Ai| and 1^2 1 remains constant 
for increasing values of 7, so that increasing 7 does not re- 
duce the uncertainty associated with the single-point distribu- 
tions of Ai,2- This indicates that synchronization is not here 
the result of reduced randomness due to the increase of time 
and length scales observed in Fig. |l]. However, the larger 
is 7 the larger becomes the mutual information, approach- 
ing its maximum possible value (/ — H{\Ai\) = 77(1^21)) 
as 7 — > 1. An additional quantitative measurement of syn- 
chronization is given by the linear correlation coefficient p = 
mi\\A2\) ~ k\AmM\)){vaT{\A^\)var{\A2\))-^l^ with 
var{x) being the variance of x. This coefficient, plotted as 
a function of 7 in Fig. I (right), is negative indicating that 
when 1 increases, |742| decreases, and vice versa. 

Our quantitative indicators of synchronization, / and p, ap- 
proach their maximum absolute values as 7 1. We also 
observe that the regime of coupled STI disappears for 7 > 1. 
In fact, the generalized synchronization observed, manifested 
by the tendency of the space-time signals towards the func- 
tional relation + |^2p = M' is probably related to the 
fact that for 7=1 this is an attracting manifold for homo- 
geneous states: Writing ^ in terms of i?^ = + 1^2 P 
and X = ai'ctan(|Ai|/|A2|), it is immediate to see that ho- 
mogeneous solutions for 7 = 1 are i?^ = /i and x arbitrary. 
To understand the preference for these solutions it is instruc- 
tive to look at the transient dynamics starting from random 
initial conditions: R{x, t) has a very fast evolution towards 
R — ^/JI with no regime of STI existing at any time. Dur- 
ing this fast evolution, the phase t) covers almost com- 
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pletely the range of its possible values. The late stages of 
the dynamics are characterized by a spatial diffusion of the 
phase xi^it) until it reaches a space-independent arbitrary 
value x{x, t) = xo- In a \Ai\ vs \A2\ dotted plot as the ones 
in Fig. ^ this is visualized by a cloud of points quickly ap- 
proaching = fi and then collapsing into a single point. 
Runs with different random initial conditions lead to different 
Xo- An underlying reason for the special dynamical behavior 
at 7 = 1 is the separation of time scales for R and x- For 
7 1, the zero wavenumber components of R and x have a 
nonzero driving force and they compete dynamically but, at 
7=1, x(fc = 0) is a marginal variable, while R{k = 0) is 
strongly driven. As a consequence, R relaxes quickly towards 
R^ ~ fi. Once R becomes space-homogeneous, the differ- 
ent wavenumber components of A12 are decoupled and the 
zerowavenumber solution wins by diffusion of x- 
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FIG. 3. Left: Entropy of | Ai | (□) and | A2 1 (A) and their mutual 
information I (o) as a function of 7. Right: Correlation coefficient p 
of I All vs 1 as a function of 7. 

In some of our simulations the STI regime has been ob- 
served to disappear for a coupling smaller than 7=1, but this 
seems to be a consequence of finite-size effects: The size of 
the laminar portions of Fig.(|l]) increases with the coupling 7. 
When this size becomes similar to system size, one of the sta- 
ble plane waves can occupy the whole system, thus prevent- 
ing any further appearance of defects and STI. For a given 
initial condition, with parameters a = 0.2 and f3 = —1.4, 
and a system size L = 512 the STI regime was seen to dis- 
appear at 7 = 0.85. As soon as the system size was doubled 
the STI regime reappeared again. By reducing system size to 
L = 256 the STI regime disappeared for smaller 7. The con- 
clusion from this an other numerical experiments is that STI 
exists for all 7 < 1 in the same range of parameters as it exists 
in the single CGLE, with time and length scales diverging as 
7 approaches 1, where STI disappears. 

In summary we have described a regime of synchronized 
STI dominated by the space-time synchronization of localized 
structures. Synchronization is measured by a mutual informa- 
tion and a correlation parameter that take their absolute maxi- 
mum value at the boundary between weak and strong coupling 
7 = 1. Beyond this boundary (7 > 1) STI disappears but the 
strong, coupled system dynamics cannot be described in terms 



of a single dominant amplitude. 
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